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FOREWORD 


The  research  work  in  this  report  was  performed  by  the 
University  of  Southampton,  Hampshire,  England,  for  the  Aero  Acoustics 
Branch,  Vehicle  Dynamics  Division,  AF  Flight  Dynamics  Laboratory, 
'Jright-Patterson  Air  Force  Base,  Ohio  under  Contract  Hr  AF  bl(052)-33S 
This  research  is  part  of  a  continuing  effort  to  obtain  data  and 
techniques  for  defining  important  vibration  *nd  acoustic  phenomena 
in  flight  vehicles  which  is  part  of  the  P.esearch  and  Technology 
Division,  Air  Force  Systems  Command's  exploratory  development 
program.  Tbs  Department  of  Defense  Pro-Tam  Element  number  is 
6.2L.05.36 A,  "Space  Flirht  Dynamics".  This  work  was  performed  under 
Project  Kr  IS^O,  "Dynamic  Problems  in  FUrht  Vehicles"  and 
Task  Nr  137005,  "Methods  of  Noise  Prediction,  Control,  and  Measurement 
Mr.  D.  I,  Smith  and  later  Mr.  ?.  H.  Hermes  of  the  AF  Flight  Dynamics 
lahoratorv  ware  the  project  engineers. 

This  work  was  administered  by  the  Office  of  Aerospace 
Research,  D3/F  through  its  European  Office  in  -(russels,  Belgium. 
Assistance  was  provided  by  Ha /or  0,  J.  Kanei,  Jr.  and 
Capt.  R.  C .  Coupland,  Jr.  of  the  ROAR, 


ABSTRACT 


Experimental  data  are  presented  on  the  r.m.s.  levels,  frequency 
spectra,  probability  distributions  and  space-time  correlations  of  the 
pressure  field  as  a  whole,  and  in  narrow  frequency  bands,  for  turbulent 
layers  on  smooth  surfaces.  General  qualitative  effect  of  a  surface 
discontinuity  on  the  r.m.s.  pressure  and  frequency  spectrum  are  given. 

An  empirical  representation  of  the  space-time  correlation  pattern  of  the 
pressure  field  is  developed  for  the  purposes  of  structural  response 
calculations.  A  general  program  for  the  study  of  a  mechanism  of 
boundary  layer  turbulence  is  presented  along  with  a  discussion  of  the 
testing  and  analyzing  instrumentation.  An  analytical  method  is  pre¬ 
sented  for  the  prediction  of  the  power  spectral  density  of  panel  res¬ 
ponse  to  boundary  layer  and  Biren  excitation.  An  empirical  equation  is 
developed  relating  the  acoustic  power  output  from  a  unit  area  of  a  (large) 
surface  subject  to  turbulent  boundary  layer  flow  to  the  aerodynamic  flow 
parameters. 


This  report  has  been  reviewed  and  is  approved. 
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1.  Introduction 


This  report  is  intended  to  summarize  the  work  carried  out  under  Contract 
No.  AF  6l(052)-358  over  the  period  March  I960  to  March  1962. 

Results  which  have  already  been  included  in  Contract  Technical  Notea 
are  briefly  summarized  with  mention  of  any  additional  results  since  obtained. 

This  applies  mainly  to  the  characteristics  of  the  turbulent  boundary  layer 
pressure  field.  Other  topics  which  have  not  been  fully  reported  on  previously, 
such  as  the  measurements  planned  for  the  laminar-turbulent  transition  region 
and  the  calculation  of  response  of  panels  to  turbulent  boundary  layer 
excitation  are  treated  In  a  little  more  detail. 

Data  have  been  obtained  on  the  r.m.s.  levels,  frequency  spectra, 
probability  distributions  and  space-time  correlations  of  the  pressure  field 
as  a  whole,  and  in  narrow  frequency  bands,  for  turbulent  boundary  layers  on 
smooth  surfaces,  and  some  effects  of  surface  discontinuities  have  also  been 
investigated.  The  work  has  reached  a  stage  where  the  statistical  character¬ 
istics  of  the  wall  pressure  field  of  a  turbulent  boundary  layer  can  be 
fairly  accurately  defined  in  terms  of  the  flow  parameters,  certainly  sufficiently 
well  to  provide  a  realistic  representation  of  the  excitation  of  a  panel  or 
structure  subjected  to  this  field.  Qnpirical  formulae  representing  the 
pressure  field  have  been  derived  for  this  purpose,  and  have  also  been  used 
to  obtain  a  figure  for  the  level  of  radiated  sound  which  would  result  from 
turbulent  flow  over  a  rigid  boundary  surface. 

The  study  of  the  wall  pressure  fluctuations  is  now  being  extended  to 

the  region  where  transition  from  laminar  to  turbulent  flow  occurs:  development 

of  the  measuring  apparatus  is  now  sufficiently  well  advanced  to  allow 

Manuscript  released  by  the  authors  in  June  .1963  for  publication  as  an 
RTD  Technical  Documentary  Report, 
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experimental  work  to  proceed  and  preliminary  results  are  included  in  this 
report. 

The  calculation  of  panel  response  to  random  loading  has  been 
programmed  for  a  digital  computer.  Estimates  of  the  response  of  a  3^"  square 
test  panel  to  turbulent  boundary  layer  excitation  have  been  made  and  similar 
calculations  are  now  being  made  to  obtain  the  response  of  the  same  panel  to 
various  forms  of  excitation  in  the  siren  tunnel.  Measurements  of  panel 
response  are  proceeding  in  parallel  with  this  work. 

Experiments  to  measure  the  levels  of  sound  transmitted  by  structural 
vibrating  resulting  from  flow  excitation  have  been  planned  to  follow  panel 
response  measurements. 

2.  Wall  Pressure  Field  Associated  with  Boundary  Laver  Flow 

2.1  Pressure  Field  of  a  Ful ly-Turbulent  Boundary  Layer 

2.1.1  Flow  Facilities  and  Instrumentation  for  Fluctuating  Pressure 
Measurements 

2.1.1. 1  Flow  Facilities 

Early  measurements  were  made  on  the  turbulent  boundary  layers  on  the 
walls  of  the  6"  x  2,V"  wind  tunnel  ^  and  the  vertical  2"  water  tunnel  ^ 

At  the  present  time  measurements  are  being  made  in  the  9"  x  6”  wind  tunnel 

II 

which  is  very  similar  in  construction  to  the  6"  x  2J-"  tunnel. 

2. 1.1.2  Pressure  Transducers 

Various  forms  of  pressure  transducer  using  piezoelectric  discs  of 
lead  zirconate-titanate  ceramic  as  sensitive  element  have  been  developed, 
with  crystal  diameters  ranging  from  0.030”  to  0.170".  The  construction  and 
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shock-tube  calibration  of  these  instruments  has  been  described  previously^ ' ^ 

Transducers  using  two  crystals  and  electrical  adding  circuits  as  a 
method  of  cancelling  unwanted  vibration  pick-up  have  been  built.  The 
per forma nee  of  one  of  these  is  referred  to  in  Section  2.2.3* 

2.1.2  Statistical  Properties  of  the  Wall  Pressure  Field 

2. 1.2.1  R.M.S.  Pressure 

Measured  values  of  r.m.s.  pressure  have  been  related  to  surface 

(2) 

friction  and  extraoolated  to  zero  transducer  size'  This  resulted  in  a 
value  of  p’/'C  D  =  2.5  at  low  speeds.  The  effects  of  increasing  flow  Mach 
number  seems  to  be  to  produce  a  slow  increase  in  the  value  of  p */tJ  to  about 
3  at  M  =  1.5.  Tests  in  the  9**  x  6"  wind  tunnel  should  provide  additional 
information  on  Mach  Number  effects. 

2. 1.2.2  Frequency  Spectrum.  Space-Time  Correlations  and  Convection 
Speeds 

The  results  of  the  experimental  measurements  of  these  quantities  have 
been  previously  presented  in  Reference  2,  where  empirical  representations  of 
the  data  were  also  given.  The  power  spectrum  of  the  pressure  fluctuations 
was  obtained  as 


which  is  consistent  with  convection  of  a  virtually  frozen  pressure  pattern 
with  a  space  correlation  in  the  flow  direction  of 


<p(x,z,t)p(x  ,zft)> 


-0.08e 


...(2.2) 


The  auto-correlation  of  the  pressure  field  referred  to  axes  moving  at  the 
convection  speed  was  given  as 
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•  •  • 


(2.3) 


Additional  experimental  work  carried  out  since  these  relationships  were 
first  obtained  has  in  general  supported  the  numerical  magnitudes  chosen.  The 
2"  water  tunnel  measurements  showed  a  rate  of  decay  of  the  pressure  field 
considerably  greater  than  was  obtained  in  other  work  and  represented  by 
Equation  (2.3).  (See  Fig.  24  of  Reference  2).  These  measurements  have  been 
examined  in  some  detail  and  corrections  applied  to  the  space-time  correlation 
values  for  the  effects  of  equipment  noise  and  a  low-level  acoustic  standing 
wave  in  the  tunnel  working  section.  The  resulting  values  are  in  very  close 
agreement  with  Equation  (2.3). 

In  Reference  2  it  was  suggested  that  in  view  of  the  lack  of  information 

the  form  of  the  space-correlation  in  the  lateral  direction  should  be  taken 

to  be  the  same  as  that  in  the  flow  direction.  Some  preliminary  measurements 

of  the  lateral  correlations  have  been  made  in  the  9"  x  6”  wind  tunnel  with 

a  flow  speed  of  about  500  ft/sec.,  but  so  far  the  minimum  transducer  separation 

it 

which  has  been  examined  is  about  5  S  •  At  this  spacing  the  correlation 
coefficient  is  very  nearly  zero  and  appears  to  be  slightly  negative.  This 
indicates  that  the  scale  of  the  field  in  the  lateral  direction  is  not 
greatly  different  to  the  longitudinal  scale. 

The  convection  speed  of  the  wall  pressure  pattern  obtained  from  water 
tunnel  measurements  increased  with  separation  of  the  measuring  points  from 
0.7Uffl  at  to  0.8511^ at  ?/  ~20.  A  similar  effect  has  now 

been  observed  in  the  9M  x  6"  wind  tunnel  where  the  convection  velocity  appears 
to  be  in  the  region  of  0.711^  for  Vs*  ~  5  increasing  to  0.80^  at  r/s*  /V  10 
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and  remaining  constant  at  this  value  for  separations  up  to  25  £*.  The  effect 
seems  to  indicate  that  the  small  scale  components  of  the  pressure  field 
result  from  rapidly  decaying  small  motions  in  the  low  mean  velocity 
region  of  the  boundary  layer  relatively  close  to  the  wall  while  the  large 
scale  components  are  associated  with  slowly  decaying  large  scale  motion  in 
the  outer  region  of  the  layer. 

2. 1.2. 3  Probability  Distribution 

An  examination  of  the  pressure  signals  from  the  water  tunnel  boundary 

(21 

layer  has  shown  that  the  amplitude  distribution  is  very  nearly  Guassian  . 

2. 1.2.4  Effects  of  Surface  Irregularities 

Results  of  measurements  of  r.m.s.  pressures  and  frequency  spectra  in 
the  region  of  step  discontinuities  in  the  wall  of  the  2"  water  tunnel  have 
been  reported  in  Reference  2.  Generalization  of  the  results  is  not  possible, 
although  some  trends  were  indicated.  For  example  the  energy  in  the  low 
frequency  region  of  the  spectrum  increased  progressively  with  the  height  of 
the  discontinuity,  but  attempts  to  obtain  a  non-dimensional  representation 
of  the  effect  in  terms  of  the  height  of  the  discontinuity  were  not  successful. 

A  probability  analysis  of  signals  from  the  disturbed  regions  showed 
that  the  relative  frequency  of  large  amplitude  excursions  tended  to  increase 
as  the  measuring  point  approached  the  discontinuity. 

Increase  of  r.m.s.  pressure  of  up  to  twice  the  smooth  wall  value  were 
observed. 

2.1.3  Empirical  Representation  of  the  Wall  Pressure  Field 

To  provide  a  self-consistent  representation  of  the  space-time 

correlation  pattern  of  pressure  field  for  purposes  of  structural  response 

calculations  the  following  empirical  representation  based  on  the  existing 

(2) 

experimental  data  was  suggested'  . 
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<p(x,z,t)p(x  +  f  ,Z  +  t  ,t  +r  »  P(f,  f,X) 


<p2> 


<Tp2> 

=  PM(  f-Uc'C  ,?,*) 


> 


pM( r. r.-o  p„(?.f »q)  pm(q.o>^) 
o2>  '  <V>  <P5> 


=  PMff 

where  Pj^f.S’t'O  is  the  space-time  correlation  of  the  pressure  field 
frame  of  reference  moving  at  the  convection  velocity  Uc«  The  form  of 
(?»  ?)  an<i  Pyt  ('•')  and  the  numerical  values  chosen  were 


<?.?>  = 


•rK2(f2+  ?2) 


.-n*2?2 


with  Kg*  =  1/4.5 


r  =  1  AL  =  -0.08 
r  =  10  A10  -  1.08 


*  t 

1 


n  =  1  BL  =  0.127 
n  =  15  B15  =  0.180 
n  -  200  B200  -  0.693 


Equations  (2.1)  to  (2.3)  above  are  based  on  these  relations. 
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2.1.4  Space-Time  Correlation  of  Filtered  Pressure  Signals 


For  calculations  of  the  response  of  a  structure  to  boundary  layer 
mode 

excitation  using  the  normal/approacf\  the  space-time  correlation  of  the 

pressure  field  in  a  narrow  frequency  band  embracing  the  particular  normal 

mode  is  required.  In  essence  this  information  is  contained  in  the  overall 

correlations  of  the  field  and  can  be  obtained  by  transformation,  but  it  can 

also  be  measured  directly.  Measurements  of  correlation  in  frequency  bands 

(4) 

were  made  first  by  Harrison  but  these  results  seemed  to  indicate  a  unique 

dependence  of  the  correlation  on  ^?/Uc*  Such  a  relation  leads  to  scale 

anomalies  at  low  frequencies.  Additional  measurements  were  made  in  the  6"  x  2jr" 

wind  tunnel ^  and  the  form  of  the  results  to  be  expected  in  the  case  of  e 

slowly  decaying  pressure  field  as  encountered  in  a  slowly  growing  boundary 

(3) 

layer  was  also  investigated  . 

It  was  concluded  in  Reference  3  that 

(1)  the  frequency  band  space-time  correlation  in  the  flow  direction 
is  periodic  in  time  delay  and  spatial  separation,  undamped  in  time,  but 

with  spatial  damping  depending  on  spatial  separation  and  convection  velocity. 

(2)  the  variation  in  correlation  amplitude  with  spatial  separation 
in  the  stream  direction  can  be  derived  directly  from  the  moving  frame 
overall  autocorrelation  of  the  pressure  field. 

A  variation  of  this  form  removes  the  scale  anomaly  at  low  frequencies 
and  on  the  basis  of  the  relations  given  in  Equations  (2.5)  and  (2.6)  it  was 
found  that  the  correlation  amplitude  could  be  expected  to  fall  from  unity  to 
0.2  with  spatial  separating  increasing  from  zero  to  about  36  $*.  The  derived 
rate  of  spatial  decay  was  found  to  be  a  tolerable  representation  of  the 
measured  frequency  band  correlations  (Fig.  7  of  Reference  3)  at  least  for 
frequencies  in  the  energy  contained  range  of  the  spectrum. 
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2.2  Wall  Pressure  Fluctuations  in  the  Laminar-Turbulent  Transition  Region 

2.2.1  Aim  of  Experiments 

The  aim  of  these  experiments  is  to  continue  the  studies  on  the  basic 
mechanism  of  boundary  layer  turbulence  by  investigating  the  properties  of  the 
turbulent  processes  found  in  the  laminar-turbulent  transition  region,  with 
particular  reference  to  wall  pressure  fluctuations. 

The  transition  region  has  been  selected  for  study  since  it  is  here 
that  the  generation  of  turbulence  first  occurs  in  a  field  with  known  properties, 
namely,  the  laminar  boundary  layer.  It  is  hoped,  therefore,  that  some  insight 
will  be  gained  into  the  orocess  of  turbulence  production,  and  hence  ultimately, 
the  mechanism  by  which  the  turbulent  boundary  layer  is  sustained. 

To  make  these  studies  use  is  being  made  of  a  hot -wire/microphone 
combination.  Selection  of  only  those  parts  of  the  flow  it  is  desired  to  study 
will  be  made  by  another  hot-wire  driving  a  specially  designed  piece  of  apparatus 
(see  below).  In  this  way,  examination  of  turbulent  spots  at  all  stages  of 
their  development  may  be  made. 

There  have  been  some  delays  in  this  programme  due  to  the  difficulty  of 
setting  up  the  instrumented  transition  region,  and  to  the  design  and  construction 
of  the  special  equipment  needed  to  study  intermittent  phenomena.  These 
difficulties  have  been  surmounted,  and  some  preliminary  results  are  presented. 

2.2.2  Experimental  Programme 

As  it  is  now  understood,  the  process  of  transition  from  laminar  to 
turbulent  flow  occurs  gradually  over  a  distance  in  the  downstream  direction 
due  to  the  superposition  of  distinct  turbulent  patches  or  spots.  These 
turbulent  spots  are  believed  to  arise  from  the  growth  of  small  disturbances 
in  the  laminar  boundary  layer  which  at  some  crucial  stage  in  their  development 
form  into  a  limited  region  of  turbulence.  This  spot  diffuses  as  it  is 
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convected  downstream  by  the  mean  flow  until  is  loses  its  separate  existence 
through  its  interaction  with  other  spots. 


Over  this  region,  where  both  laminar  and  turbulent  flow  regimes  are 
to  be  found,  the  degree  of  "spottiness"  is  described  by  the  intermittency 
of  the  flow.  The  intermittency  may  be  defined  either  as  the  fraction  of  time 
for  which  the  flow  is  laminar  (the  "laminar  intermittency",  denoted  here  by 
G),  or  the  fraction  of  time  for  which  the  flow  is  turbulent  (the  "turbulent 
intermittency",  denoted  by  Y  ),  As  the  latter  is  the  one  in  general  use  in 
the  literature,  the  term  intermittency  will  hereafter  be  used  for  turbulent 
intermittency. 

It  is  intended  to  study  the  turbulence  in  the  spots  in  the  transition 
region.  In  particular,  it  is  hoped  to  study  the  pressure  fluctuations 
associated  with  the  transition  region,  the  character  of  the  turbulence  in 
the  spots,  and  the  position  in  the  layer  where  the  turbulence  ap  iears  first. 
This  last  aspect  may  be  of  importance  in  relation  to  the  observed  convection 
of  the  pressure  field  of  a  turbulent  layer  at  0.3  of  the  free  stream 
velocity,  since  it  has  been  suggested  that  this  effect  is  due  to  phenomena 
near  the  surface  excited  in  conditions  in  the  outer  part  of  the  boundary  layer. 

Recent  advances  in  knowledge  of  transition  have  been  made  by 
Schubauer,  Klebanoff  and  others ^  using  a  vibrating  ribbon  to  artificially 
induce  transition.  It  is  not  planned  at  the  moment  to  use  this  technique, 
but  to  study  the  statistical  characteristics  of  the  flow  with  natural 
transition. 

The  basic  configuration  for  the  instrumentation  is  shown  in  Fig.  11. 
Studies  of  the  field  will  be  made  using  a  microphone  and  hot-wire  (marked 
"active  wire"),  and  the  outputs  from  each  will  be  analyzed  both  separately 
and  also  as  a  joint  record  by  correlation  techniques.  In  order  that 
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the  turbulence  in  the  region  of  the  microphone  may  be  limited  to  newly 
developed  spots,  another  hot-wire  placed  upstream  of  the  microphone  (marked 
"monitor  wire"  will  be  used  to  operate  a  device  which  will  shut  off  the 
microphone  signal  if  the  turbulence  is  passing  over  the  monitor  wire.  Although 
only  one  channel  is  treated  in  this  way,  the  effect  after  correlation  is 
as  though  both  active  hot-wire  and  microphone  outputs  were  processed  by  the 
monitor  wire. 

2.2.3  Instrumentation  of  the  Transition  Region 

In  order  to  make  a  study  of  the  transition  region  with  acoustic 
instrumentation  it  is  necessary  to  have  some  control  over  the  position  at  which 
transition  occurs,  since  the  exact  position  is  very  sensitive  to  external 
flow  conditions  which  vary  from  day  to  day,  and  it  is  not  in  general  Dossible 
to  move  the  microphones  about  readily.  Basically,  in  a  wind  tunnel,  the  option 
is  between  a  boundary  layer  on  a  surface  specially  mounted  for  that  purpose, 
or  using  the  tunnel  wall. 

At  the  time  these  measurements  were  first  envisaged,  the  only  low 
noise  tunnel  immediately  available  was  an  induced  flow  one  with  a  working 
section  measuring  6"  x  2.5"  which  had  been  in  use  for  turbulent  boundary 
layer  measurements.  The  speed  in  the  working  section  of  this  tunnel  was  fixed 
at  about  500  f.p.s.  by  a  downstream  sonic  choke.  As  the  flow  had  been 
stabilized  by  placing  boundary  layer  trips  at  the  entrance  it  was  not  possible 
to  study  transition  on  the  working  section  walls,  and  so  a  plate  was  mounted 
in  the  working  section  for  this  purpose.  The  high  frequencies  characteristic 
of  the  flow,  and  the  very  limited  space  available,  made  the  use  of  the  piezo¬ 
electric  crystal  transducers  previously  used  in  this  tunnel  for  pressure 
measurements  imperative.  It  was  not  possible  in  the  space  to  incorporate 
any  form  of  mechanical  vibration  isolation  for  the  crystal,  so  a  vibration 


cancelling  technique  was  evolved  using  two  crystals.  The  two  crystals  were 
mounted  rigidly  on  the  plate  and  one  of  them  was  screened  from  the  pressure 
field;  their  electrical  outputs  were  combined  to  cancel  the  acceleration 
effects.  Unfortunately,  the  presence  of  the  plate  excited  strong  acoustic 
resonances  in  the  tunnel,  and  as  it  would  have  required  extensive  modification 
to  the  tunnel  to  avoid  them  this  arrangement  was  abandoned.  The  method  of 
vibration  compensation  of  the  pressure  transducers  is  of  some  interest,  however, 
and  the  degree  of  success  can  be  judged  from  Fig.  1,  where  the  suppression 
of  the  plate  modes  Is  evident. 

At  this  stage  a  new  tunnel  of  a  similar  type  to  the  one  mentioned 
above,  but  with  a  working  section  measuring  9"  x  6",  was  complete.  It  was 
thought  that  a  suitable  transition  region  could  be  established  on  the  wall 
of  this  tunnel  at  the  beginning  of  the  working  section,  the  position  of 
transition  being  controlled  by  sucking  away  various  amounts  of  the  boundary 
layer.  However,  the  flow  tended  to  become  turbulent  upstream  of  the  suction 
slots,  so  that  satisfactory  control  of  transition  could  not  be  obtained.  It 
was  not  considered  feasible  with  the  conditions  present  to  re-stabilize  the 
boundary  layer  once  it  had  become  turbulent. 

It  was  decided,  therefore,  to  revert  to  study  of  the  transition  region 
on  a  flat  plate,  and  to  adapt  a  low  speed  tunnel  specially  for  the  work, 
replacing  the  existing  entrance  with  a  new  inlet,  working  section  and 
diffuser,  and  altering  the  fan  gearing.  The  existing  aluminium  sheet  diffuser 
was  treated  with  "Aquaplas",  and  as  much  of  the  new  structure  as  possible 
was  made  with  acoustically  inert  fibre  board.  The  resulting  tunnel  offers 
a  speed  of  up  to  100  f.p.s.  in  the  working  section.  A  diagram  of  the 
arrangement  is  shown  in  Fig.  2,  and  of  the  working  section  in  Fig.  3. 

The  working  surface  is  a  brass  plate,  highly  polished,  with  a  flat 
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upper  surface  and  a  6°  wedge  cut  on  the  leading  edge.  A  slight  waviness 
is  visible  to  the  eye,  but  the  surface  is  smooth  to  the  touch.  A  small 
hole,  diameter  0.030",  is  drilled  7"  from  the  leading  edge,  giving  access  to 
a  cavity  for  a  condenser  microphone  (BrQel  and  Kjaer  type  4134)  which  is 
mounted  on  the  underside  of  the  plate,  and  recessed  to  within  0.029"  of 
the  top  surface.  Control  of  boundary  layer  transition  is  effected  by 
adjusting  the  pressure  gradient  in  the  working  section  by  means  of  flexible 
wall  liners. 


2.2.4  Turbulence  Monitoring  and  Gating  Equipment 

Fundamental  to  the  concept  of  these  experiments  is  the  ability  of  the 
instrumentation  to  distinguish  between  signals  characteristic  of  laminar  and 
turbulent  flow.  The  way  in  which  the  apparatus  for  doing  this  function  is 
described  below. 

The  core  of  this  apparatus  is  a  detecting  (or  rectifying)  system  to 
which  a  variable  bias  may  be  applied.  To  illustrate  the  working,  suppose 
for  simplicity  that  the  signal  to  be  examined  is  a  sine  wave,  as  shown  in 
Fig.  4  and  we  wish  to  know  for  what  fraction  of  the  neriod  T  the  amplitude 
is  above  some  arbitrary  level,  y.  It  will  be  seen  that  at  time  t^  the 
amplitude  becomes  greater  than  y,  and  stays  so  until  time  t2S  the  appropriate 
fraction  is 


since 


G  =  (tj-tjJA  -*[l 

tj  =  sin-1y 

t2  *  IT  -  sin-1y 


...(2.7) 


The  output  of  an  ideal  device  for  a  given  setting,  y,  would  be  given 
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by  Equation  2.7,  and  varies  from  G  =  0  (level  never  above  y  for  y  >  1 )  to 
G  =  0.5  (half  wave  rectification)  for  y  =  0,  when  the  apparatus  is  examining 
positive  parts  of  the  signal.  In  practice  it  is  necessary  that  an  increasing 
(decreasing)  signal  exceed  (fall  below)  the  setting  y  by  some  small  amount 
before  the  apparatus  can  be  made  aware  of  it.  Errors  introduced  by  this 
effect  are  greatest  at  the  peak.  This  and  other  similar  effects  due  to 
finite  operating  times  can,  however,  be  reduced  by  careful  design.  Some 
idea  of  the  magnitude  of  possible  errors  can  be  obtained  from  Table  1. 


TABLE  1 

Error  due  to  Non-Ideal  Operation  of  Trigger  Circuit 
for  Positive  Going  Portion  of  Sine  Waves 


Necessary  excess  in/ 
fraction  of  signal 

0.2 

n 

n 

0.8 

0.95 

0.01 

0.003 

0.004 

0.005 

0.009 

0.018 

0.05 

0.007 

0.008 

0.010 

0.018 

0.035 

0.10 

0.033 

0.038 

0.049 

0.088 

0.175 

Error  given  in  fraction  of  whole  period. 


If  this  basic  device  is  presented  with  an  intermittent  signal  -*nd 
the  level  y  made  just  high  enough  not  to  be  reached  by  the  instrument 
noise  (see  Fig.  6),  then  the  device  will  "fire"  at  the  peaks,  but  there  will 
be  a  dead  zone  around  the  axis,  equal  to  2y,  in  which  it  will  shut  off.  To 
overcome  this,  it  is  necessary  to  build  in  a  "memory"  which  holds  the  device 
on  for  long  enough  to  ensure  that  there  is  no  further  critical  peak  to  follow. 
The  last  super-critical  peak  will  be  followed  then  by  a  fallacious  extension 
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to  the  time  of  the  spot.  This  process  is  shown  in  Fig.  6, where  the  unmodified 
gating  signal  and  signal  with  time  delay  are  illustrated  for  a  hypothetical 
input  signal. 

The  circuit  diagram  is  shown  in  ^ig.  8,  and  block  diagram  for  the 
circuit  in  Fig.  9.  The  signal  input  (marked  "Monitor  input")  is  amplified 
and  passes  to  a  signal  inverter,  which  has  a  two-channel  output,  one  being 
the  same  as  the  input  signal  and  the  other  the  input  signal  inverted.  This 
is  done  so  that  both  positive  and  negative  portions  of  the  signal  can  be 
examined,  the  detecting  circuit  being  sensitive  to  positive  going  signals 
only.  After  further  amplification  each  part  of  the  signal  passes  to  the 
detector  system,  which  is, in  essence,  a  diode  to  which  a  variable  bias  may 
be  applied,  thus  passing  the  proportion  of  the  signal  above  the  critical 
(bias)  level.  The  resulting  clipped  signal  is  amplified  so  that  the  input 
to  the  Schmitt  Trigger  has  sharp  rise  characteristics.  The  Schmitt  Trigger 
produces  square  pulses  as  wide  as  the  clipped  peaks  remaining  from  the 
detector.  These  pulses  pass  to  the  Time  Delay,  which  imparts  an  exponential 
decay  to  their  trailing  edges,  the  rate  of  decay  being  variable  by  simply 
altering  the  appropriate  capacitances.  The  result  is  to  effectively  widen 
the  pulses  which  are  presented  to  the  second  Schmitt  Trigger,  the  output  of 
which  is  used  to  operate  a  gating  circuit. 

Any  signal  may  be  applied  for  gating  ("Active  input"),  including  the 
monitor  signal  itself.  As  arranged,  the  gate  is  shut  when  the  critical 
levels  are  exceeded  in  the  detector,  and  hence  is  normally  open.  Thus,  the 
result  of  passing  any  intermittent  signal  through  the  gating  system  will  be 
to  remove  those  parts  of  it  which  occur  at  times  when  the  monitor  signal 
exceeds  the  critical  level. 

A  counting  system  has  also  been  provided  for  use  with  the  circuit 
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for  quantitative  measurements  of  intermittency .  In  this  case,  the  input  to 
the  gate  is  a  200  kcps  sine  wave.  The  output  is  passed  through  three 
binary  counters,  which  reduce  the  effective  frequency  to  12.5  kcps,  before 
operating  a  timed  decade  counter.  In  this  way  the  fraction  of  the  time  for 
which  the  monitor  signal  is  sub-critical  is  simply  the  ratio  of  the  counts 
with  and  without  monitor  input. 

The  last  mentioned  facility  provides  an  accurate  means  of  assessing 
the  accuracy  of  the  device,  since  by  varying  the  clip  level  in  the  detector 
it  is  possible  to  determine  how  much  cf  the  given  signal  is  greater  than  a 
given  amplitude  -  effectively  a  probability  analyser.  Thus  for  a  sine  wave, 
the  fraction  is  given  by  Equation  2.7,  while  a  similar  evaluation  may  be  made 
for  white  noise.  The  results  of  actual  tests  are  shown  in  Figs.  5  and  7.  It  is 
seen  that  the  operation  of  the  equipment  is  satisfactory. 

2.2.5  Results  Obtained  to  Date 

As  this  equipment  was  only  made  operational  a  short  time  before  this 
writing  it  has  not  been  possible  to  accumulate  any  extensive  experimental 
results.  However,  some  preliminary  results  nave  been  taken,  and  are 
presented  below. 

2.2.5. 1  Calibration 

The  initial  calibration  of  the  microphone  showed  a  flat  response 
to  about  2000  c.p.s.  after  which  acoustic  resonances  produced  a  rather  peaky 
and  generally  rising  response  to  about  10,000  c.p.s.  The  application  of 
damping  material  in  the  form  of  steel  wool  in  the  probe  bore  has  reduced  the 
peaks,  although  the  gradual  rise  is  still  present.  All  the  same  the  response 
curve  is  flat  to  within  -  4  db  up  to  10  kcps.,  which  is  above  the  expected 
cut-off  frequency  for  the  spectrum  of  the  turbulence  given  by  the  Strouhal 
number  relation 
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leading  to 


f  S*/Uoo~  0,1 

£^6000  c.p.s.  for  =  100  f.p.s. 

S*  =  0.020" 

The  background  noise,  which  is  dominated  by  discrete  fan  frequencies 
at  about  150  c.p.s.  and  (to  a  lesser  extent)  at  300  c.p.s.,  is  at  an 
acceptable  level  of  about  94  db,  although  the  low  frequencies  do  show  up  in 
the  correlation  coefficient  obtained  from  hot-wire  and  microphone  signals, 
and  will  be  eliminated  by  filtering  when  necessary.  The  free  stream  turbulence 
level,  as  measured  with  the  new  low-noise  transistor  hot-wire  gear,  is  about 
0.3*6  of  the  free  stream  mean  velocity,  and  as  this  level  does  not  seem  to  be 
having  any  adverse  effect  on  the  transition  phenomenon,  no  attempt  is  being 
made  to  reduce  it  at  present. 

To  examine  different  parts  of  the  transition  region,  it  is  necessary 
to  alter  the  pressure  gradient  slightly  with  the  tunnel  liners.  While  this  gives 
the  microphone  effectively  a  variable  x  co-ordinate  relative  to  the  beginning 
of  the  transition  region,  it  also  has  some  effect  on  the  length  of  the  actual 
zone  over  which  transition  occurs.  An  estimate  of  the  extent  of  the  transition 
region  is  to  be  had  if  it  is  assumed  that  the  plot  of  intermittency  v.  distance 
along  the  plate  resembles  an  integrated  Gaussian  Error  function.  There  is 
experimental  support  for  this  assumption.  The  intermittency  measured  at 
two  points  4.2"  apart  was  0.138  and  0.778.  For  this  case  turbulence  can  be 
expected  for  1^  of  the  time  4.5"  upstream  of  the  leading  station,  and  90&  at  3.5" 
downstream  of  the  trailing  station,  although  the  assumption  of  Gaussian  form 
is  almost  certainly  in  error  at  this  distance.  However,  the  extent  of  the 
region  is  probably  of  the  order  of  about  9". 
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2. 2. 5. 2  Intermittency  Measurements 


Fig.  12  shows  intermittency  measured  at  various  stations  in  the  boundary 
layer  Tor  two  stations  in  the  transition  region,  one  where  the  layer  is  mainly 
laminar  and  the  other  where  the  layer  is  mainly  turbulent.  Though  it  is 
inevitable  with  such  measurements  that  there  should  be  rel=tively  large 
scatter,  both  curves  are  quite  well  defined.  It  will  be  seen  that  there 
exists  a  marked  variation  in  the  measured  intermittency  across  the  boundary 
layer,  maximum  values  being  observed  at  about  0.15  S  and  0.65  S  ,  where  S  is 
the  thickness  of  the  undisturbed  laminar  layer.  If  one  oart  of  the  layer 
predominates  in  the  production  of  turbulent  spots,  and  it  is  assumed  that 
the  spots  grow  uniformly  as  they  are  convectod  downstream,  it  may  bo  shown 
that  the  observed  intermittency  in  this  region  will  be  higher  than  elsewhere. 

It  is  possible,  therefore,  that  one  or  both  of  the  maxima  observed  in  Fig.  12 
correspond  to  areas  of  turbulence  production. 

2. 2. 5. 3  Pressure  Measurements 

Measurements  have  also  been  made  of  the  pressure  fluctuations,  using 
a  hot-wire  placed  immediately  over  the  microphone  to  evaluate  the  intermi ‘'tency. 
The  hot-wire  signal  is  to  be  preferred  for  this  purpose  to  the  microphone 
output  as  the  signal/noise  ratio  is  much  better,  but  it  has  the  disadvantage 
that  (as  can  be  seen  from  Fig.  12  )  the  intermittency  measured  ■'t  a  fixed 
distance  from  the  surface  will  only  approximately  represent  the  value  of  the 
wall  pressure  intermittency;  the  value  obtained  in  this  way,  however,  can  be 
expected  to  be  fairly  representative  of  the  pressure  intermittency.  The 
microphone  output  was  integrated  over  a  period  of  30  seconds  and  divided  by 
this  intermittency  to  arrive  at  a  value  for  the  r.m.s.  wall  oressure 
associated  with  the  turbulent  spots. 

Measurements  with  full  turbulence  at  a  flow  speed  of  70  feet  per 
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second  gave  a  pressure  level  of  104  db  re  0.0002  dyne/cm  .  For  this  measure¬ 
ment  the  microphone  was  calibrated  against  a  standard  noise  source.  This 
makes  the  ratio 

p'/q  =  9.5  x  10~3 

where  p'  is  the  r.m.s.  value  of  the  pressure  fluctuations,  and  q  is  the 
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dynamic  pressure,  which  is  the  same  as  that  obtained  by  Harrison  ' .  This 
preliminary  figure  is  subject  to  revision,  but  indicates  no  significant  difference 
from  the  fully  developed  turbulent  boundary  layer. 

The  pressure  associated  with  the  spots  for  different  intermittencies 
(different  effective  downstream  stations)  is  shown  in  Fig.  10,  together  with 
the  accuracy  limits.  As  the  intermittency  approaches  zero,  (fully  laminar), 
both  the  integrated  value  of  the  pressure  signal  and  the  divisor  (  Y  )  become 
very  small.  If  the  error  in  the  r.m.s.  pressure  P  is  p,  and  in  the  measured 
intermittency  Y  is  g,  then  the  quotient  is 


Since  approximately  {>  =  kP 

then  p-  ■  P  1  ;  kp/y  f  (1  t  *  T  3  ) 

>  l!*  '»  r  7 

p  and  g  are  fixed  in  any  given  situation  so  that  the  limit  lines  vary 
hyperbol ical ly  with  ^  .  The  limits  are  given  for  2%  error  in  P,  -,nd  3%  in 
V  .  Both  are  rather  optimistic  estimates,  but  it  is  seen  that  the  resulting 
lines  nearly  enclose  the  pressure  points  from  fully  turbulent  down  to  about 
3C£i  turbulent.  Agreement  is  poorer  at  very  low  values  where  the  error  is 
of  the  order  of  the  mean  value.  The  mean  value  used  in  Fig.  10  may  be  made 
to  coincide  with  the  fully  turbulent  value  (  Y  =  1.0)  by  using  a  different 
value  of  intermittencyj  as  discussed  in  Section  2. 2.5.2  above,  the  precise 


intermittency  which  is  applicable  is  still  in  doubt.  Fig.  11  then  indicates 
that  the  pressure  field  associated  with  the  turbulent  spots  is  similar  to 
that  of  fully  developed  turbulence. 

2. 2, 5. 4  Continuation 

Work  is  proceeding  for  the  time  being  with  the  single  wire  and 
microphone  combination  discussed  above.  In  particular,  further  work  is 
being  done  on  the  variation  of  intermittency  across  the  boundary  layer,  and 
on  the  variation  of  amplitude  distribution,  which  has  been  observed  to  be 
markedly  skew.  It  is  proposed  not  to  proceed  with  two  wire  work  until  the 
possibilities  of  the  present  configuration  have  been  exploited;  the 
necessary  facilities  for  doing  this  are  already  available. 

3,  Theoretical  Response  of  Panel  to  Turbulent  Boundary  Layer  Excitation 

3.1  Re suits  to  Date 

The  work  done  so  far  has  been  based  on  the  results  derived  by  Powell ^ 
for  the  response  of  structures  to  random  excitations,  and  solutions  are, 
therefore,  obtained  as  the  sum  of  solutions  in  each  ot  a  multitude  of  normal 
modes.  The  limitation  of  this  approach  is  that  response  in  each  mode  has 
to  be  derived  separately,  with  the  consequence  that,  in  general,  a  large  amount 
of  computation  is  needed  to  get  even  an  approximate  evaluation  of  the  panel 
displacements  in  quite  a  narrow  frequency  range. 

As  yet  the  range  of  freauencies  covered  is  limited  to  the  region 
0  -  1000  c/s.  which  includes  the  first  15  panel  modes. 
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3.2  Basis-  of  Calculations 


3.2. 1  Equations  used 


Equation  (13a)  of  Reference  6  for  the  power  spectral  density  of  panel 
displacement  can  be  written  in  the  form 


W^f)  _  y, 

.  r 


r2(K)AI  j2  (f) 


?  /  hr|  i  vs| 

=  —  I  J  ~f(<  ,  )fT(s)l/rs(s')  d’s  d's< 

AI  .A  JA  ...(3.2 

=  J  l^r2(s)  d2j>  j  A  =  Panel  Area  ...(3.3 


AI  j2Nrs(f) 


... (3.1) 


where  n  1 

.r  ^  /  r  \  _ 


j  Nrs(f)  =  a^Tfl  AI 


...(3.2) 


Panel  Area 


...(3.3) 


Here  (f ) ,  ,‘tt'p(f )  are  power  spectral  densities  of  oanel  displacement 
z  at  jl,  and  exciting  pressure  p  respectively,  at  a  frequency  f  c/s. 

*»  p(f;i,s' j  -r)  is  the  cross-power  spectral  density  of  the  excitation, 

ifr r(s)  is  the  mode  shape  considered  evaluated  at  point  js. 

pf  (  S  »  s  *  T  )  =  P—  (^>^»T)  is  normalised  cross-power  spectral  density 


of  the  excitation 


“>  -  2f  f 


K  -  9sl 


Y  =  / Y  |  e  1  is  panel  impedance  for  mode  r. 
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3.2.2  Panel  Impedance 


The  form  used  here  is 


Mr  f 


(wT2  -  w  2)2  +  ft-2  W  4' 


■-*’1  f'-(f)2] 


£ 


...(3.4) 


Here  ,  a  "hysteretic"  damping  coefficient,  is  taken  as  0.005,  but  the 
actual  value  taken  has  little  significance,  except  for  frequencies  near  to 
the  natural  frequency  for  the  mode  considered,  provided  only  that  it  is  of 
the  right  order  of  magnitude. 

3.2.3  Effect  of  Cross  Terms 

As  stated  in  the  Section  3.1  we  are  confining  our  activities  for  tho 
moment  to  the  low  frequency  range,  so  that  it  seems  to  be  justifiable  to 
assume  that 


i •> 

will  not  be  too  close  to  zero,  and,  hence,  that  P^ ( £  »  C  >T  )  will  be 
appreciably  less  than  P^(J,^,0).  This  fact,  combined  with  the  fact  that 
the  integral  of  a  oroduct  of  modal  functions  over  part  of  the  panel  area 
will  usually  be  less  than  the  integral  of  the  square  of  either  of  the  modal 
functions  over  the  same  area,  leads  us  to  assume  that 


r  /  s 

The  validity  of  this  assumption  for  the  particular  correlation 


...(3.6) 
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function  assured  is  one  of  the  matters  that  is  to  be  investigated  later. 
Equation  (3.1)  is  thus  used  in  the  reduced  form 


4  ^/(R)  AI  J2Nrr(f)  ...(3.7) 

Vf)  W2 


3.3  Correlation  Function 
3.3.1  General  Form 

In  Reference  3  the  following  form  for  »f  »T  )  for  a  slowly 

decaying  pressure  pattern  is  deriveds 


pf  ($,$,*)  -  nmr? 


it-o 


"■4;  •  °) 


C0SW(T"  *_) 

C  C 


8) 


where  P„  (r)  is  the  autocorrelation  function  referred  to  axes  moving  with 

MT 

convection  velocity  Uc»  whilst  the  term 


NW1*  (  uc  *\) 
nmr^  (u;  » °) 

is  taken  to  be  the  same  as  the  overall 
So  (3.8)  is  taken  as 


lateral 


space  correlation 


(0,  9  ). 


VM’t>"pM!f(0’5)p Mt(^)  C08‘*^-  u;) 

(See  Reference  3,  Section  4  from  which  this  is  taken). 

i.e.  Pf(|,^,0)  =  x  PL0N3*  ...(3.10) 

PLAT’  PLONG  are  overa11  lateral  space  correlation  and  longitudinal  space 
correlation  in  required  narrow  frequency  bands  respectively. 
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3.3.2  Assumed  Form  for  Boundary  Laver  Tunnel  Correlation  Function 
Because  of  its  analytical  simplicity,  the  form  which  has  been  assumed 
initially  for  P  (r  )  is  that  taken  by  Dyer  in  Reference  7.  However,  the 
computer  programme  can  be  extended  to  include  correlation  functions  such  as 
that  proposed  in  Reference  2. 


where 

Here 


P*,T  <  *  ) 


,-M  /• 


U 


00 


30 

500  ft/sec 


5* 


0.1" 


so  that  6  =  0.  5  m.sec. 
Thus  we  can  take 


...(3.11) 


P  =  P 

Lore  mt 


=  exp 


f  S  ^  cos  “5 

Uc 

l'  ^  il cos  4  t) 


...(3.12) 


Jq  =  separation  in  longitudinal  direction  for  first  zero  of  P, 
-  log  decrement  based  on  first  (negative)  peak. 


Then 


L0N3 


i.e. 

for 

=  2  | 

io 

2tf 

* 

2f  \± 

% 

Uc 

L 

V 

— 

100 

ft. 

9  0 

f 

= 

1200 

ins 

...(3.13) 
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2  So 


Uc® 


S  L  =  1000 


...(3.14) 


There  Is  very  little  experimental  information  on  the  value  of  It  will 

be  assumed  that 


$  =  2.5  i  *  =  0.25  ins 

o 

and  the  maximum  negative  value  is  0.05. 

So  we  take 


PLAT  =  6XP 


l‘ix  iJ 
2  ?ol 


cos  2t 

i — 

|  *. 

. 

V,  i 

=  3.0 


... (3.15) 

...  (3.16) 
...(3.17) 


where 

Equations  (3.12)  -  (3.17)  inclusive  substituted  into  (3.10)  give  the  assumed 
form  for  the  correlation  function. 

3.3.3  Assumed  Form  for  the  Correlation  Function  _for_  the  Siren  Tunnel 

Here  we  assume  P.  =  1  ...(3.18) 

LAI 

(Reference  8  Figure  14), 

and  assuming  that  for  a  given  frequency  the  effective  part  of  the  pressure 
field  is  a  cosine  wave  propagating  at  the  soeed  of  sound,  we  get 

...(3.19) 

is.  ...(3. 20 ) 


*L  =  0 


£  =  3300  ins 

>  o  f 


so  that 


Pf  (  ^  >  o)  *  plqn3 

with  values  given  in  Equations  (3.19)  and  (3.20). 


...  (o..l) 
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3.4  Evaluation  of  j2Nrp(f) 

9 

In  Reference  9  a  relationship  is  derived  for  j  (f )  as  defined  in 
Equation  (3.2)*,  for  a  correlation  function  of  the  form 


'aL^  1 

e  cos  bjj 


...(3.22) 


with  ^Js)  taking  the  form 

'f‘I  (l)  5^fmn^  =  sin  mirx  sin  n,ry 


...(3.23) 


The  mode  shape  assumed  here  is  appropriate  to  having  the  panel  edges  simply 

supported  but  there  are  indications  that  the  effect  of  the  plate  boundary 

2 

conditions  on  the  integrated  value  of  j  jgrr(f)  Is  small  for  the  fundamental 
mode,  and  gets  progressively  smaller  as  the  order  of  the  mode  increases,  so 
that  the  error  introduced  by  the  assumption  is  probably  small. 

Using  relations  similar  to  Equations  (3.22)  and  (3.23),  Reference 
9  derives 


(f) 


1 

(mn-r2)2 


...(3.24) 


where 


0  = 


sin  bL 


...(3,25) 

...(3.26) 


*This  definition,  rather  than  that  used  by  Powell,  is  used  to  eliminate  the 
effects  of  any  arbitrary  coefficients  arising  in  the  choice  of Tfr(s). 
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...(3.27) 


...(3.28) 


... (3.29) 


where  all  symbols  except  numerical  constants  are  understood  to  have  suffix 
and  where  i i  -  L,  T. 

m  =  m,  mj  =  n. 


The  correspondence  between  Equation 
is  given  by 


3.5  Evaluation  of  J^pJ^ 


(3.22)  and  Equation  (3.12)  and  (3.16) 


1 

...  (3.30) 

o 

CM 

S  T 

...  (3.31) 

2<o 

3.5.1 

From  Equation  (3.5) 

I *rl  ■  TV*'!  -(3-32) 

where  Mp  =  ^  d*J,  =  ml  ...(3.33) 

m  =  mass/unit  area  of  panel. 

Thus  the  evaluation  of  j Yr |  2  centres  on  the  evaluation  of  the  natural 
frequencies  of  the  panel  in  each  of  the  modes  considered. 
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3.5.2  Evaluation  of  tal 

_  r 

The  natural  frequency  of  the  panel  is  derived  using  the  work  reported 
in  Reference  10.  The  actual  boundary  conditions  (i.e.  fixed  or  supported 
depending  on  the  panel)  are  used  here,  as  distinct  from  the  assumption  made 
for  Equation  (3.23). 

3.6  Final  Equation 

We  now  have  all  that  is  required  for  the  evaluation  ofT&z(f)  since 
for  Equation  (3.8) 


tbz(f) 

Vf) 


m.n 


32N«) 


...(3.34) 


mn 


where  !mn  =  M/^mn^  ^  ...(3.35) 

'A 

and  the  ^^(s),  « m  are  those  dictated  by  the  actual  plate  boundary  condition. 
3.7  Panel  Results 

The  panel  being  evaluated  at  the  moment  is  an  0.005"  thick  steel 
panel,  3.5"  square,  and  the  panel  is  fastened  to  a  solid  member  in  such  a 
way  that  it  is  reasonable  to  take  it  as  having  fully  fixed  edges. 

The  process  indicated  above  has  been  carried  through  for  this  oanel, 
the  lowest  15  modes  having  been  considered  for  the  response  in  the  boundary 
layer  tunnel,  but  only,  so  far,  the  lowest  mode  for  the  siren  tunnel 
case. 


The  results  of  the  computations  are  shown  in  Fig.  13. 
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4.  Response  of  Panels  to  Boundary  Laver  Excitation  (Experimental) 

4.1 

The  study  of  structural  response  to  turbulent  boundary  layer  excitation 
is  of  interest  in  relation  to  two  major  engineering  problems}  one  is  the 
problem  of  noise  in  and  around  aircraft  and  submarines  where  a  major 
contribution  results  from  the  action  of  the  structure  as  a  transducer  in 
converting  the  hydrodynamic  pressures  in  the  boundary  layer  into  acoustic 
pressures  and  the  other  is  that  prolonged  excitation  by  the  boundary  layer 
may  cause  structural  failure  due  to  fatigue.  In  an  investigation  into 
either  of  these  problems  it  is  desirable  to  have  a  knowledge  of  the  response 
of  simpler  elements  of  a  structure,  such  as  panels,  to  random  excitation. 

Panel  response  to  boundary  layer  excitation  has  been  studied 
theoretically  by  several  authors  using  either  panel  normal  modes  or  the 
concept  of  coincidence  (i.e.  frequency  and  wavelength  matching  between 
exciting  force  and  panel  response).  However,  up  to  date  there  has  been 
little  or  no  experimental  verification  of  the  accuracy  of  estimation  provided 
by  either  of  these  approaches  and  hence  it  is  the  object  of  the  current 
investigation  to  obtain  data,  by  measuring  the  response  of  panels  to  boundary 
layer  excitation  in  air,  which  might  provide  (a)  justification  for  the 
calculation  procedures  and  (b)  a  clearer  insight  into  tha  mechanism  of  the 
response  of  a  structure  to  this  type  of  excitation.  The  problem  is  complex 
and  the  model  structures  studied  initially  must  of  necessity  be  simple  in 
form  so  that  they  are  amenable  to  calculation.  For  this  reason  the 
measurements  are  being  made  on  panels  of  square  or  rectangular  plan-form, 
having  fixed  edges  and  having  no  longitudinal  or  transverse  stiffeners. 

In  order  to  interpret  the  measurements  fully,  the  vibrational 
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characteristics  of  the  specimens  must  be  known.  For  this  reason  the 
frequencies  and  shapes  of  the  normal  modes  of  the  panels  are  determined 
experimentally  using  discrete  frequency  excitation.  To  complete  the 
investigation,  the  response  of  the  panels  under  test  is  being  calculated 
theoretically,  as  described  in  Section  (3)  above,  for  comparison  with  the 
experimental  values. 

4.2  Preliminary  Measurements  in  the  6"  x  2^-”  Wind  Tunnel 

The  initial  work  on  panel  response  to  boundary  layer  noise  excitation 
was  carried  out  in  the  supersonic  working  section  of  the  6”  x  2^-"  tunnel. 

It  was  intended  to  investigate  possible  coincidence  effects  at  supersonic 
convection  velocities  of  the  boundary  layer  pressure  field,  with  particular 
reference  to  the  generation  of  Mach  waves  by  supersonic  flexural  ripples  in 
the  panel.  This  study  proved  impracticable  in  this  small  scale  flow  and 
when  Mr.  D.J.M.  Williams  left  the  group,  this  work  was  discontinued  in 
favour  of  the  panel  displacement  measurements  in  the  9"  x  6"  tunnel. 

The  work  did,  however,  highlight  one  of  the  major  difficulties  in 
these  investigations,  namely  the  construction  of  the  panel.  When  using  thin 
sn vets  of  metal  it  is  extremely  difficult  to  obtain  a  panel  which  is  flat 
and  yet  under  zero  or  negligible  tension.  The  panel  was  formed  from  a  sheet 
of  brass,  O.OCft"  thick,  attached  to  a  6"  diameter  plug,  the  shape  of  the 
vibrating  panel  being  determined  by  a  3^-"  square  hole  cut  through  the  centre 
of  the  plug.  The  first  panel  was  stuck  to  the  plug  with  epoxy  resin  but 
the  adhesive  failed  during  trial  runs  in  the  tunnel.  Later  panels  were 
attached  to  the  plug  by  clamping  strips  parallel  to  the  edges  of  the  square 
hole.  This  method  presented  difficulties  in  fixing  the  brass  sheet  at  the 
corners  of  the  cut-out  and  introduced  unknown  and  probably  non-uniform 
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tension  in  the  panel. 

4.3  Test  in  the  9"  x  6"  Wind  Tunnel 

The  measurements  in  the  9"  x  6"  tunnel  are  concerned  with  the  dis¬ 
placements  of  panels  subjected  to  boundary  layer  excitation  and  hence  are 
of  a  rather  different  nature  to  those  discussed  above.  The  panels  being 
tested  are  similar  to  those  used  in  the  6"  x  2^"  tunnel  in  that  they  are 
formed  from  a  thin  sheet  of  metal  attached  to  a  6"  diameter  plug  with  the 
hole  in  the  plug  again  determining  the  shape  of  the  vibrating  region.  The 
panels  currently  being  used  are  3^-"  square  and  either  0.005"  or  0.010"  thick. 
The  material  used  has,  however,  been  changed  from  brass  to  steel,  so  that 
the  panels  can  be  excited  magnetically  for  determination  of  their  normal 
modes  of  vibration. 

Because  of  the  possibility  of  introducing  unknown  tensions  in  the 
panels,  and  hence  creating  unpredictable  effectr.  in  their  response ^  the 
clamping  procedure  adopted  in  the  earlier  tests  was  not  continued.  Instead, 
recourse  was  had  to  epoxy  resins  for  bonding  the  steel  sheets  to  the  plugs. 

It  was  found  that,  with  care,  the  adhesives  could  withstand  the  static 
pressure  differentials  imposed  on  them  when  placed  in  the  wind-tunnel  wall. 

To  obtain  a  good  bond  the  surfaces  are  completely  degreased  by  immersion  in 
a  lC/o  solution  of  sodium  metasilicate  for  a  period  of  10  minutes.  The 
brightness  of  the  surfaces  is  then  enhanced  by  immersing  in  an  85&  phosphoric 
acid  solution  for  2  minutes. 

Whilst  the  epoxy  resin  is  setting,  the  plug  and  steel  sheet  are  held 
together  in  a  jig  so  that  the  panel  remains  flat.  This  clamping,  which 
was  not  used  when  assembling  the  first  0.005"  panel,  does  not  cause  panel 
tensioning  and  the  jig  is  removed  when  the  resin  has  hardened.  To  avoid 
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differential  thermal  expansion  between  the  plug  and  panel  due  to  localised 
heating,  the  resin  is  cured  under  normal  room  temperature.  This  does 
however  have  the  disadvantage  of  increasing  the  curing  time  and  slightly 
decreasing  the  strength  of  the  bond. 

The  construction  of  a  typical  oanel  is  shown  in  Fig.  14.  The  two 
clamping  strips,  along  sections  of  the  plug  perimeter,  are  positioned  after 
the  resin  has  set  and  hence  do  not  tension  the  panel.  Their  purpose  is  purely 
precautionary  and  is  to  prevent  the  panel  from  being  carried  into  the  tunnel 
in  the  event  of  a  failure  of  the  adhesive  during  tunnel  operation.  The 
region  around  the  clamping  strips  is  built-up  with  resin  to  oresont  a  surface 
which  is  flush  with  the  tunnel  wall. 

4.4  Normal  Modes 

The  normal  modes  of  the  panels  are  determined  from  panel  response  to 
discrete  frequency  excitation.  The  natural  frequencies  are  obtained  from 
measurements  of  amplitude  and  phase  (relative  to  the  exciting  force)  of  the 
panel  vibrations  and  the  mode  shapes  are  observed  visually  when  sand  particles 
are  spread  ovc-r  the  oanel  surface. 

The  frequencies  of  the  low  order  modes  were  investigated,  using 
electromagnetic  excitation,  for  the  0.CO5"  panel  made  first.  The  values 
obtained  were  found  to  be  considerably  higher  than  those  estimated 
theoretically  (Status  Report  No.  8).  In  particular,  the  measured  funda¬ 
mental  frequency  was  in  the  range  160  c.p.s.  to  260  c.p.s.  (subject  to 
positioning  of  equioment  as  discussed  below)  whilst  the  calculated  value  was 
141  c.p.s. 

It  is  believed  that  the  discrepancy  is  due  to  the  fact  that  the  panel, 

when  assembled,  is  not  flat  but  has  an  appreciable  curvature.  This  bow  is 

very  difficult  to  control  but  a  satisfactory  degree  of  flatness  can  probably 
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be  achieved  by  the  modified  method  of  construction  described  in  Section  4.3. 
Some  approximate  calculations  based  on  comparisons  between  the  natural 
frequencies  of  spherical  caps  as  derived  by  Reissner  in  Reference  11 
and  those  for  flat  circular  plates  (Reference  12)  indicate  that  the  degree 
of  bowing  to  be  expected  can  have  an  effect  on  the  panel  of  similar  magnitude 
to  that  measured  above.  As  the  natural  frequencies  of  flat  square  plates 
of  side  a,  and  flat  circular  plates  of  diameter  a  differ  by  less  than 
3%,  it  is  reasonable  here  to  take  the  results  for  circular  plates  and  apply 
them  directly  to  square  ones. 

Closer  examination  of  the  measurements  showed  that  the  "apparent 
natural  frequencies"  depended  markedly  on  the  relative  positions  of  the 
exciter,  panel  and  proximity  transducer  used  to  measure  the  vibrations.  This 
was  most  critical  when  exciter  and  transducer  were  both  at  the  panel  centre 
and,  for  these  conditions,  the  observed  variation  of  panel  "fundamental" 
frequency  with  distance  of  transducer  from  panel  is  shown  in  Fig.  15. 
Unfortunately,  oefore  this  phenomenon  could  be  fully  explored,  the  panel 
was  damaged  and  had  to  be  replaced.  The  investigation  will,  however,  continue 
on  the  second  c.anel  using  electromagnetic  and  acoustic  excitation  -  the 
latter  being  of  use  only  for  modes  for  which  the  generalised  force  is  non¬ 
zero. 

4.5  Boundary  Laver  Excitation 

The  panel  is  mounted  in  the  side  wall  of  the  9"  x  6"  Induced  flow 
wind  tunnel.  The  static  pressure  in  the  tunnel  working  section  is  below 
atmosoheric  and  in  the  sub-sonic  section  the  pressure  differential  is  2.3 
to  2.4  p.s.i.,  dependong  on  position}  this  has  a  marked  effect  on  the 
vibration  of  the  thin  panels  placed  in  the  tunnel  wall. 

The  displacement  power  spectrum  for  the  first  0.005"  panel  under  a 
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pressure  differential  of  2.34  p.s.i.  was  shown  in  Fig. 2  of  Status  Report 
No.  8.  The  main  purpose  of  these  measurements  was,  however,  a  comparison 
of  the  amplitudes  of  the  panel  vibration  and  the  background  vibration  of  the 
traverse  gear  and  wind  tunnel.  It  was  found  that  the  latter  vibration  was 
at  least  20  db  below  the  former  throughout  the  spectrum  and  would  therefore 
not  influence  measurements  of  panel  vibration. 

The  outer  face  of  the  panel  has  now  been  enclosed  by  a  sealed  box  in 
which  the  static  pressure  can  be  reduced  to  that  of  the  tunnel.  The  background 
vibration  is  till  of  the  same  order  of  magnitude  as  previously  measured  but 
the  overall  panel  vibration  amplitude  has  Increased  by  a  factor  of  about  4. 
Panel  displacement  spectra  have  been  measured  at  several  positions  on  the 
surface  of  a  second  0.005"  panel  and  results  for  three  of  these  positions 
are  shown  in  Fig.  16.  The  pressure  differential  across  the  panel  during  the 
vibration  measurements  was  0.003  p.s.i. 

The  spectra  show  the  presence  of  several  panel  normal  modes  and  are 
in  this  respect  considerably  different  in  character  from  those  previously 
measured  under  a  greater  pressure  differential.  Fig.  16  clearly  shows  the 
variation  of  the  relative  magnitude  of  the  response  in  each  mode  with 
variation  of  position  on  the  surface  (i.e.  as  the  vibration  pick-up  chnnges 
position  in  relation  to  the  modes  and  antimodes  of  the  mode),  for  example, 
the  peak  at  790  c.p.s.  (Fig.  16  (a))  is  barely  discernible  in  Fig.  16(b)  but 
has  returned  again  in  Fig.  16(c). 

The  natural  frequencies  of  this  panel  have  not  yet  been  determined 
cXpcriiTcnta  lly  and  hence  the  particular  modes  present  in  the  above  spectra 
cannot  yet  be  identified.  The  frequency  of  the  predominant  peak  (295  c.p.s.) 
does  not  however  correspond  with  the  calculated  fundamental  frequency 
(141  c.p.s. )  of  the  panel. 
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5.  Sound  Radiation  from  a  Turbulent  Boundary  Laver  on  a  Rigid  Surface 

On  the  basis  of  the  correlation  pattern  given  by  Equations  (2.4)  to 

(2.6)  the  acoustic  power  output  from  unit  area  of  a  (large)  surface  subject 

(2 ) 

to  turbulent  boundary  layer  flow  was  found  to  be  given  by'  ' 


W 


<  f»u»6  s/»»3j 


...(5.1) 
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OPERATION  OF  MONITORING 
EQUIPMENT  ON  SINE  WAVE 
INPUT. 
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FIGURE  3:  CIRCUIT  JIAGUAi.  Faft  TUBBULEKT  iXCiKUBIBG  EQUIPMENT 
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PI6URE  »1.  BASIC  CONFIGURATION  FOR  MEASUREMENTS. 


LiSURE  12.  VARIATION  OF  INTERHITTEHCV  ACROSS  SOUHftWY  LAVER. 
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FjCURE  If  EFFECT  OF  DISTANCE  OF  PROXIMITY  PICK-UP  FROM  SURFACE  ON 
NATURAL  FREQUENCY  OF  PANEL  EXCITEO  ELECTRO  MAGNETICALLY. 
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